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Abstract 


Ships  appear  in  synthetic  aperture  radar  (SAR)  ocean  imagery  as  bright  targets  against  the  ocean 
clutter  background.  Envisat  Advanced  SAR  (ASAR)  Alternating  Polarization  (AP)  mode  data 
acquisitions  over  the  Strait  of  Dover  and  the  Strait  of  Gibraltar,  two  high-density  shipping  regions 
with  Automatic  Identification  System  (AIS)  coverage  via  AISLive,  have  provided  a  large 
database  of  validated  ship  signatures.  For  each  validated  ship  signature,  several  metrics  were 
computed  including  the  clutter  statistics,  the  total  radar  cross  section  of  each  ship,  and  the  ship 
signature  length.  These  metrics  were  evaluated  in  terms  of  the  polarization,  the  local  incidence 
angle,  the  validated  ship  length,  the  target  to  clutter  ratio,  and  the  ship  aspect  angle.  This  is  a  very 
rich  data  set  that  has  yielded  new  observations  and  insights  to  ship  detectability  by  SAR.  These 
results  could  have  bearing  on  the  design  of  future  SAR  modes  or  SAR  missions,  and  could  be 
used  to  improve  the  performance  of  ship  detection  software. 


Resume 


Sur  les  images  de  l’ocean  acquises  au  radar  a  synthese  d’ouverture  (RSO),  les  navires 
apparaissent  comme  des  cibles  brillantes  sur  Tarriere-plan  de  fouillis  d’echos.  Des  donnees 
acquises  au  moyen  du  RSO  evolue  (ASAR)  de  l’Envisat  en  mode  de  polarisation  altemee  (PA) 
sur  le  Pas  de  Calais  et  le  detroit  de  Gibraltar,  deux  regions  a  grande  densite  d’activite  maritime 
couvertes  par  le  Systeme  d’identification  automatique  (S1A  en  direct,  par  TAlSLive),  ont  foumi 
une  grosse  base  de  donnees  de  signatures  de  navires  validees.  Pour  chaque  signature  de  navire 
validee,  plusieurs  mesures  ont  ete  calculees  dont  les  statistiques  sur  le  fouillis  d’echos,  la  section 
efficace  en  radar  de  chaque  navire  et  la  longueur  de  la  signature  du  navire.  Ces  mesures  ont  ete 
evaluees  en  termes  de  polarisation,  d’angle  d’incidence  local,  de  longueur  validee  du  navire,  du 
rapport  cible  sur  fouillis  et  d’angle  d’aspect  du  navire.  II  s’agit  d’un  tres  riche  jeu  de  donnees  qui 
a  fourni  de  nouvelles  observations  et  perspectives  sur  la  detectabilite  des  navires  au  moyen  du 
RSO.  Ces  resultats  pourraient  etre  importants  pour  la  conception  de  tuturs  modes  ou  missions 
RSO  et  pourraient  etre  utilises  pour  T  amelioration  du  rendement  des  logiciels  de  detection  de 
navires. 
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Executive  summary 


Validation  of  ship  signatures  in  Envisat  ASAR  AP  mode  data 
using  AISLive:  Data  acquisition,  processing,  and  analysis 
results 

Vachon,  Paris  W.;  Wolfe,  John;  DRDC  Ottawa  TM  2008-005;  Defence  R&D 
Canada  -  Ottawa;  March  2008. 

Introduction:  Ships  appear  in  synthetic  aperture  radar  (SAR)  ocean  imagery  as  bright  targets 
against  the  ocean  clutter  background.  Envisat  Advanced  SAR  (ASAR)  Alternating  Polarization 
(AP)  mode  data  acquisitions  over  the  Strait  of  Dover  and  the  Strait  of  Gibraltar,  two  high-density 
shipping  regions  with  Automatic  Identification  System  (AIS)  coverage  via  AISLive,  a 
commercial  AIS  provider,  have  provided  a  large  database  of  validated  ship  signatures  that  can  be 
used  to  address  the  polarization  dependence  of  ship  detectability.  The  AIS  data  were  in  turn 
validated  through  the  Internet  Ships  Register  (ISR). 

The  analysis  methods  are  similar  to  those  used  to  analyze  previously  acquired  RADARSAT-1 
Fine  and  ScanSAR  Narrow  B  mode  data  sets  of  ships  validated  by  AISLive  data. 

Results:  For  each  validated  ship,  several  metrics  were  computed  using  the  VUSAR  image 
analysis  tool.  These  included  the  clutter  statistics,  the  total  radar  cross  section  (RCS)  of  each 
ship,  and  the  ship  signature  length.  These  metrics  were  evaluated  in  terms  the  polarization,  the 
local  incidence  angle,  the  validated  ship  length,  the  target  to  clutter  ratio,  and  the  ship  aspect 
angle.  The  following  is  apparent: 

•  The  cross-polarization  ocean  clutter  is  at  the  Envisat  ASAR  noise  floor,  so  any  target 
contrast  metrics  would  be  relative  to  the  instrument  noise  floor  rather  than  the  ocean  clutter; 

•  The  background  ocean  clutter  may  be  characterized  as  K-distributed; 

•  For  the  ship  targets,  the  HH  and  VV  total  RCS  are  comparable  while  the  cross-polarization 
total  RCS  is  about  10  dB  smaller  than  the  co-polarization  total  RCS;  models  for  the  total 
RCS  as  a  function  of  ship  length  are  provided; 

•  For  co-polarization,  the  ship  target  contrast  increases  with  increasing  incidence  angle,  while 
for  cross-polarization  the  contrast  is  more-or-less  independent  of  incidence  angle;  the  co¬ 
polarization  ship  contrast  is  larger  for  cross-polarization  if  the  incidence  angle  is  smaller 
than  33°; 

•  The  total  RCS  variability  was  successfully  modelled  by  fitting  the  normalized  total  RCS  to 
several  well  known  probability  density  functions;  the  RCS  variability  is  smaller  for  cross¬ 
polarization  ship  signatures; 

•  The  ship  signature  length  agreed  with  the  validated  ship  length  for  60.1%  of  the  co¬ 
polarization  ships  and  71.5%  of  the  cross-polarization  ships;  for  co-polarization,  the  rate  of 
underestimation  of  ship  length  is  higher;  for  cross-polarization,  the  rate  of  overestimation  of 
ship  length  is  higher;  there  tends  to  be  more  overestimates  of  ship  length  for  azimuth 
travelling  ships; 
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•  Ship  wake  signatures  appear  in  the  cross-polarization  channel  as  a  ship-like  signature  that 
has  not  been  subject  to  azimuth  shifting,  as  occurs  for  the  ship  signature  itself;  these 
signatures  appear  frequently  and  may  be  due  to  scattering  from  waves  generated  off  the  bow 
of  the  ship. 

Significance:  This  is  a  very  rich  data  set  that  has  yielded  new  observations  and  insights  to  ship 
detectability  by  SAR.  These  results  could  have  bearing  on  the  design  of  future  SAR  modes  or 
SAR  missions,  and  could  be  used  to  improve  the  performance  of  ship  detection  software. 

The  use  of  AIS  data  for  ship  signature  validation  is  now  a  well-proven  methodology  that  should 
be  used  to  develop  ship  detection  models  for  other  missions  such  as  TerraS AR-X  and  Cosmos- 
Skymed.  Focussing  on  high  density  shipping  regions  would  provide  a  significant  set  of  validated 
ship  targets  for  each  image  acquired. 


IV 
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Sommaire 


Validation  of  ship  signatures  in  Envisat  ASAR  AP  mode  data 
using  AISLive:  Data  acquisition,  processing,  and  analysis 
results 

Vachon,  Paris  W.;  Wolfe,  John;  DRDC  Ottawa  TM  2008-005;  R  &  D  pour  la 
defense  Canada  -  Ottawa;  Mars  2008. 

Introduction  ou  contexte:  Sur  les  images  de  1’ ocean  acquises  au  radar  a  synthese  d’ouverture 
(RSO),  les  navires  apparaissent  comme  des  cibles  brillantes  sur  l’arriere-plan  de  fouillis  d’echos. 
Des  donnees  acquises  au  moyen  du  RSO  evolue  (ASAR)  de  l’Envisat  en  mode  de  polarisation 
altemee  (PA)  sur  le  Pas  de  Calais  et  le  detroit  de  Gibraltar,  deux  regions  a  grande  densite 
d’activite  maritime  couvertes  par  le  Systeme  d’identification  automatique,  SIA  en  direct,  par 
l’AISLive,  un  foumisseur  commercial  du  SIA,  ont  foumi  une  grosse  base  de  donnees  de 
signatures  de  navires  validees  qui  peut  etre  appliquee  a  la  determination  de  la  dependance  de  la 
detectabilite  des  navires  a  l’endroit  de  la  polarisation.  Les  donnees  du  SIA  ont  ete  a  leur  tour 
validees  d’apres  V Internet  Ships  Register  (ISR). 

Les  methodes  d’analyse  sont  similaires  a  celles  appliquees  pour  l’analyse  de  jeux  de  donnees 
RADARSAT-1  anterieurement  acquises  en  mode  fin  et  en  mode  ScanSAR  etroit  B  sur  des 
navires  et  validees  d’apres  des  donnees  du  AIS  Live. 

Resultats:  Pour  chaque  signature  de  navire  validee,  plusieurs  mesures  ont  ete  calculees  au  moyen 
de  l’outil  VUSAR  d’analyse  d’images,  dont  les  statistiques  sur  le  fouillis  d’echos,  la  section 
efficace  en  radar  (RCS)  de  chaque  navire  et  la  longueur  de  la  signature  du  navire.  Ces  mesures 
ont  ete  evaluees  en  termes  de  polarisation,  d’angle  d’incidence  local,  de  longueur  validee  du 
navire,  du  rapport  cible  sur  fouillis  et  d’angle  d’aspect  du  navire.  On  a  releve  ce  qui  suit : 

•  le  fouillis  d’echos  de  l’ocean  en  polarisation  orthogonale  est  au  plancher  de  bruit  de  l’ASAR 
de  l’Envisat  de  sorte  que  toute  mesure  du  contraste  de  cible  serait  etablie  par  rapport  au 
plancher  de  bruit  de  l’instrument  plutot  que  par  rapport  au  fouillis  d’echos  de  l’ocean; 

•  le  fouillis  d’echos  d’arriere-plan  de  l’ocean  peut  etre  caracterise  comme  presentant  une 
distribution  K; 

•  pour  les  cibles  que  sont  les  navires,  les  RCS  totales  en  HH  et  en  VV  sont  comparables  alors 
que  la  RCS  totale  en  polarisation  orthogonale  est  d’environ  10  dB  inferieure  a  la  RCS  totale 
en  copolarisation;  des  modeles  de  la  RCS  totale  en  fonction  de  la  longueur  du  navire  sont 
foumis; 

•  en  copolarisation,  le  contraste  des  navires  cibles  augmente  lorsque  augmente  Tangle 
d’incidence  alors  qu’en  polarisation  orthogonale  le  contraste  est  plus  ou  moins  independant 
de  Tangle  d’incidence;  le  contraste  des  navires  cibles  en  copolarisation  est  plus  grand  en 
polarisation  orthogonale  si  Tangle  d’incidence  est  inferieur  a  33°; 

•  la  variabilite  de  la  RCS  totale  a  ete  modelisee  avec  succes  en  ajustant  la  RCS  totale 
normalisee  a  plusieurs  fonctions  de  distribution  de  probabilites  bien  connues;  la  variabilite 
de  la  RCS  est  moindre  pour  les  signatures  de  navires  en  polarisation  orthogonale; 
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•  la  longueur  de  la  signature  de  navire  correspondait  avec  la  longueur  validee  des  navires  a 
60,1  %  en  copolarisation  et  a  71,5  %  en  polarisation  orthogonale;  en  copolarisation,  le  taux 
de  sous-estimation  de  la  longueur  des  navires  est  plus  eleve;  en  polarisation  orthogonale,  le 
taux  de  surestimation  de  la  longueur  des  navires  est  plus  eleve;  il  y  a  tendance  a  davantage 
de  surestimation  de  la  longueur  des  navires  pour  les  navires  en  deplacement  en  azimut; 

•  sur  le  canal  de  polarisation  orthogonale  les  signatures  de  sillages  de  navires  s’apparente  aux 
signatures  de  navires  n’ayant  pas  subi  de  decalage  en  azimut,  comme  c’Est  le  cas  pour  la 
signature  du  navire  elle-meme;  ces  signatures  se  manifestent  frequemment  et  peuvent  etre 
attribuables  a  la  diffusion  engendree  par  les  vagues  soulevees  a  l’etrave  du  navire. 

Importance:  Le  jeu  de  donnees  decrit  est  tres  riche  et  a  foumi  de  nouvelles  observations  et 
perspectives  sur  la  detectabilite  des  navires  par  RSO.  Ces  resultats  pourraient  etre  importants  pour 
la  conception  de  futurs  modes  ou  missions  RSO  et  pourraient  etre  utilises  pour  L  amelioration  du 
rendement  des  logiciels  de  detection  de  navires. 

L’utilisation  des  donnees  du  S1A  pour  la  validation  des  signatures  de  navires  est  une  methode 
eprouvee  qui  devrait  etre  utilisee  pour  1’ elaboration  de  modeles  de  detection  de  navires  dans  le 
cadre  d’autres  missions  comme  les  TerraSAR-X  et  Cosmos-Skymed.  Le  fait  de  concentrer  les 
travaux  sur  les  a  grande  densite  d’activite  maritime  foumirait  un  important  jeu  de  navires  cibles 
valides  pour  chaque  image  acquise. 


VI 
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1  Introduction 


Ships  appear  in  synthetic  aperture  radar  (SAR)  ocean  imagery  as  bright  targets  against  the  ocean 
clutter  background.  Ship  signature  characterization  in  RADARSAT-1  (R-l)  and  Envisat 
Advanced  SAR  (ASAR)  imagery  can  improve  automatic  ship  detection  algorithms,  identify 
signature  metrics  for  ship  classification,  and  provide  information  for  predicting  the  ship  detection 
performance  of  future  SAR  systems. 

Recently,  R-l  (C-band,  HH  polarization)  Fine  mode  and  ScanSAR  Narrow  B  (SCNB)  mode  data 
sets  were  successfully  acquired,  along  with  AISLive  fwww.  aislive.  com)  data  that  were  used  for 
ship  signature  validation  [1],  [2],  [3].  This  document  describes  the  extension  of  the  R-l  activity 
to  include  other  polarizations  by  using  the  Envisat  ASAR  (C-band)  Alternating  Polarization  (AP) 
mode.  This  work  provides  new  insight  to  the  impact  of  polarization  choices  on  ship  detection. 

The  available  data  sets  and  a  high-level  summary  of  the  analysis  results  from  all  trials  have 
already  been  described  elsewhere  [4], 

Envisat  ASAR  AP  mode  data  acquisition  over  the  Strait  of  Dover  and  the  Strait  of  Gibraltar,  two 
high-density  shipping  regions  with  AISLive  coverage,  commenced  in  mid-November  2006  and 
finished  by  mid-February  2007.  A  total  of  42  alternating  polarization  (AP)  mode  images 
representing  all  seven  ASAR  beams  (i.e.,  IS1  through  IS7,  represented  the  smallest  to  largest 
incidence  angles)  and  including  all  three  available  dual  polarization  configurations  (i.e.,  F1F1/F1V, 
VY/VF1,  and  F1F1/VV)  were  ordered,  along  with  42  time-coincident  pairs  of  AISLive  snapshot 
data  sets  bracketing  each  ASAR  image  acquisition  time. 

The  Envisat  ASAR  images  were  obtained  through  European  Space  Agency  (ESA)  Announcement 
of  Opportunity  (AO)  project  #255  (P.W.  Vachon,  Principle  Investigator).  The  images  were 
obtained  from  ESA  as  radiometrically  calibrated,  precision  16-bit  products.  Some  of  the  products 
were  obtained  via  the  Envisat  Rolling  Archive  within  a  few  days  of  acquisition.  Others  were 
obtained  via  CD-ROM. 

Of  the  42  images  ordered,  35  ASAR  images  were  eventually  received.  The  remaining  7  images 
were  not  acquired  due  to  programming  conflicts  or  technical  issues  with  the  satellite  operation. 
All  of  the  corresponding  AISLive  data  sets  were  received  via  email  within  a  few  days  of  data 
acquisition. 

The  ships  present  in  the  35  acquisitions,  as  reported  by  AISLive,  were  validated  through  the 
Internet  Ships  Register  (ISR)  (www.ships-register.com).  The  validation  procedure  consists  of 
matching  the  AISLive  ship  position  information  with  the  closest  SAR  target  signature,  and 
performing  a  handshake  via  the  International  Maritime  Organization  (IMO)  number  to  access  the 
definitive  ISR  database  ship  specifications  (i.e.,  the  ship  length). 

Processing  of  the  AISLive  data  was  carried  out  as  described  in  [3]  to  establish  the  AlS-projected 
ship  position  at  the  Envisat  pass  time,  the  AlS-predicted  ship  position  that  accounts  for  azimuth 
shifting  of  the  AlS-projected  position,  the  actual  ship  signature  position,  the  ship  radar  cross 
section  (RCS),  and  other  ship  signature  metrics.  The  analysis  methods  are  similar  to  those  used 
to  analyze  the  previously  acquired  R-l  Fine  and  SCNB  data  sets. 
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Section  2  contains  a  summary  of  the  ASAR  and  validation  data  sets;  Section  3  describes  the 
properties  of  the  background  clutter  against  which  the  ship  signatures  appear;  Section  4  addresses 
the  ship  RCS  and  related  contrast  measures;  Section  5  models  the  variability  of  the  observed  ship 
RCS;  Section  6  compares  the  ship  signature  length  with  the  validated  ship  length;  and.  Section  7 
contains  some  examples  of  ship  wakes  observed  in  the  cross-polarization  channel  data. 
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2  Data  summary 


The  ASAR  data  were  ordered  using  Earthnet  OnLine  Interactive  (EOLI),  a  GUI-based  tool  that 
provides  access  to  the  online  ESA  catalogues  of  EO  products.  AISLive  data  sets  were  ordered 
and  received  via  e-mail  from  Lloyd’s  Register-F airplay  once  the  ASAR  orders  were  confirmed. 

Table  1  summarizes  the  ASAR/AISLive  acquisition  parameters.  A  red  row  indicates  an  ASAR 
programming  conflict  while  a  magenta  row  indicates  that  beam  IS5  was  inoperable  and  that  no 
ASAR  data  for  that  beam  mode  and  polarization  combination  were  acquired.  Each  data  set 
consists  of  one  image  and  two  AISLive  snapshots,  one  acquired  before  the  ASAR  pass  time  and 
one  acquired  after  the  ASAR  pass  time,  collected  nominally  3  minutes  apart. 

The  basic  statistics  for  the  number  of  validated  ships  in  the  4  linear  polarizations  are  shown  in 
Table  2.  The  various  parameters  in  the  table  are  described  elsewhere  [3].  The  total  number  of 
ships  available  for  each  polarization  state  has  been  highlighted  in  yellow  and  represents  high 
confidence  matches  between  the  AlS-predicted  location  and  the  target  signature  location. 

The  AISLive-reported  ship  length  was  verified  against  the  ISR  database  ship  length  via  a 
scatterplot,  as  shown  in  Figure  1.  There  are  several  different  measures  of  ship  length  (e.g.,  length 
overall,  waterline  length,  scantling  length).  Therefore,  an  error  was  defined  in  the  verification 
process  as  a  difference  between  the  AISLive-reported  length  and  the  ISR  database  length  overall 
of  30%  or  more.  The  error  rate  is  close  to  10.5%,  which  is  slightly  larger  than,  but  generally 
consistent  with,  the  error  rate  seen  in  the  previous  AISLive  data  sets.  Ship  length  errors  include 
the  ship  length  reported  in  feet  rather  than  meters,  ship  width  or  some  other  measure  appearing  in 
the  ship  length  field,  and  most  frequently,  a  reported  ship  length  of  zero.  In  a  few  cases,  a  ship 
length  was  not  available  from  the  ISR  database. 

Histograms  of  the  separation  distance  between  the  AlS-predicted  position  and  the  target  signature 
position  are  shown  in  Figure  2,  with  one  histogram  presented  for  each  polarization.  The  mean  (~ 
150  m)  and  standard  deviation  (~  75  m),  as  shown  in  the  title  of  each  sub-plot,  are  similar  for 
each  polarization. 

Histograms  of  the  ship  lengths  extracted  from  the  ISR  Database  for  the  validated  ships  in  each 
linear  polarization  are  shown  in  Figure  3.  The  mean  (~  150  m)  and  standard  deviation  (~  60  m) 
are  similar  for  each  polarization.  Note  that  the  mean  distance  between  the  AlS-prediction 
position  and  the  target  signature  position  is  about  the  same  as  the  mean  ship  size  the  database. 
This  reflects  good  performance  in  the  use  of  AIS  for  ship  signature  validation. 

Histograms  of  the  estimated  total  RCS  for  the  ships  in  each  polarization  are  shown  in  Figure  4. 
As  with  the  R-l  data  considered  previously,  there  is  a  very  wide  range  of  radar  cross  sections 
observed,  with  the  cross-polarization  RCS  being  systematically  smaller  than  the  co-polarization 
RCS. 

Validated  ships  are  not  identically  visible  in  each  available  channel  of  AP  mode  data.  This  is 
especially  true  for  smaller  incidence  angles.  Examples  of  recognizable  VH  ship  signatures  versus 
undetectable  VV  ship  signatures  at  near  incidence  angles  are  shown  in  Figure  5.  Similarly,  range- 
ward  transects  of  normalized  radar  cross  section  at  small  and  large  incidence  angles  are  shown  in 
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Figure  6.  Note  how  the  target  contrast  becomes  larger  for  larger  incidence  angles,  especially  for 
co-polarization  at  low  incidence  angle,  which  is  dominated  by  clutter. 


Table  1:  ASAR/AISLive  acquisition  parameters. 


Site 

Swath 

Pol 

Orbit 

Track 

Pass 

Date 

ASAR 

StartTime 

ASAR 
Stop  Time 

ASAR 

CentreTime 

AIS 

Snapl 

AIS 
Snap  2 

Dover 

ISl 

HH/HV 

24761 

158 

A 

24-Nov-06 

21:27:36 

21:28:51 

21:28:13 

21:27:13 

21:30:13 

VV/VH 

24990 

387 

A 

10-Dec-06 

21:24:45 

21:25:01 

21:24:53 

21:23:53 

21:26:55 

Dover 

IS2 

HH/HV 

25885 

280 

D 

11 -Feb-07 

10:22:24 

10:22:39 

10:22:31 

10:21:31 

10:24:31 

VV/VH 

25534 

430 

A 

17-Jan-07 

21:30:25 

21:30:40 

21:30:32 

21:29:33 

21:32:33 

HH/VV 

25384 

280 

D 

7- Jan-07 

10:22:24 

10:22:40 

10:22:32 

10:21:33 

10:24:33 

Dover 

IS3 

HH/HV 

24804 

201 

A 

27-Nov-06 

21:33:18 

21:33:33 

21:33:25 

21:32:26 

21:35:25 

VV/VH 

25805 

201 

A 

5 -Feb-07 

21:33:16 

21:33:31 

21:33:23 

21:32:23 

21:35:24 

HH/VV 

25305 

201 

A 

1 -Jan-07 

21:33:17 

21:33:32 

21:33:24 

21:32:24 

21:35:24 

Dover 

IS4 

HH/HV 

25849 

244 

A 

8-Feb-07 

21:39:00 

21:39:15 

21:39:07 

21:38:08 

21:41:10 

VV/VH 

25348 

244 

A 

4- Jan-07 

21:38:59 

21:39:14 

21:39:07 

21:38:06 

21:41:06 

HH/VV 

25076 

473 

A 

16-Dec-06 

21:36:09 

21:36:24 

21:36:16 

21:35:16 

21:38:16 

Dover 

1S5 

HH/HV 

25620 

15 

A 

23-Jan-07 

21:41:51 

21:42:06 

21:41:59 

VV/VH 

25255 

151 

D 

29-Dec-06 

10:05:19 

10:05:34 

10:05:26 

HH/VV 

25119 

15 

A 

19-Dec-06 

21:41:51 

21:42:06 

21:41:59 

Dover 

IS6 

HH/HV 

25162 

58 

A 

22-Dec-06 

21:47:33 

21:47:48 

21:47:41 

21:46:41 

21:49:41 

HH/VV 

24661 

58 

A 

17-Nov-06 

21:47:33 

21:47:48 

21:47:41 

21:46:12 

21:49:12 

Dover 

IS7 

HH/HV 

25434 

330 

A 

1 0-Jan-07 

21:50:24 

21:50:39 

21:50:32 

21:49:32 

21:52:32 

VV/VH 

24933 

330 

A 

6-Dec-06 

21:50:24 

21:50:40 

21:50:32 

21:49:32 

21:52:32 

HH/VV 

25935 

330 

A 

14-Feb-07 

21:50:25 

21:50:40 

21:50:32 

21:49:32 

21:52:32 

Gibraltar 

ISl 

HH/HV 

24969 

366 

D 

9-Dec-06 

10:38:07 

10:38:22 

10:38:15 

10:37:16 

10:40:15 

VV/VH 

25334 

230 

A 

3 -Jan-07 

22:06:30 

22:06:45 

22:06:37 

22:05:37 

22:08:37 

HH/VV 

25198 

94 

D 

25-Dec-06 

10:35:16 

10:35:32 

10:35:24 

10:34:25 

10:37:25 

Gibraltar 

IS2 

HH/HV 

24926 

323 

D 

6-Dec-06 

10:32:26 

10:32:41 

10:32:34 

10:31:36 

10:34:34 

VV/VH 

25105 

1 

A 

1 8-Dec-06 

22:09:20 

22:09:35 

22:09:28 

22:08:28 

22:11:28 

HH/VV 

25427 

323 

D 

10-Jan-07 

10:32:26 

10:32:41 

10:32:34 

10:31:34 

10:34:34 

Gibraltar 

IS3 

HH/HV 

25155 

51 

D 

22-Dec-06 

10:29:35 

10:29:51 

10:29:43 

10:28:45 

10:31:43 

VV/VH 

25656 

51 

D 

26- Jan-07 

10:29:36 

10:29:51 

10:29:43 

10:28:43 

10:31:43 

HH/VV 

25148 

44 

A 

21 -Dec-06 

22:15:01 

22:15:16 

22:15:08 

22:14:18 

22:17:09 

Gibraltar 

IS4 

HH/HV 

24919 

316 

A 

5-Dec-06 

22:17:52 

22:18:07 

22:17:59 

22:17:02 

22:19:59 

HH/VV 

25384 

280 

D 

7- Jan-07 

10:26:45 

10:27:00 

10:26:52 

10:25:53 

10:28:53 

Gibraltar 

IS5 

HH/HV 

25191 

87 

A 

24-Dec-06 

22:20:43 

22:20:58 

22:20:50 

22:19:50 

22:22:50 

VV/VH 

25112 

8 

D 

19-Dec-06 

10:23:54 

10:24:09 

10:24:01 

10:23:02 

10:26:02 

HH/VV 

25692 

87 

A 

28-Jan-07 

22:20:43 

22:20:58 

22:20:50 

Gibraltar 

IS6 

HH/HV 

24962 

359 

A 

8-Dec-06 

22:23:33 

22:23:48 

22:23:40 

22:22:40 

22:25:41 

VV/VH 

25341 

237 

D 

4- Jan-07 

10:21:04 

10:21:19 

10:21:11 

10:20:11 

10:23:11 

HH/VV 

25842 

237 

D 

8 -Feb-07 

10:21:04 

10:21:19 

10:21:11 

10:20:11 

10:23:07 

Gibraltar 

IS7 

HH/HV 

25069 

466 

D 

16-Dec-06 

10:18:13 

10:18:28 

10:18:20 

10:17:22 

10:20:20 

VV/VH 

25735 

130 

A 

31 -Jan-07 

22:26:24 

22:26:39 

22:26:32 

22:25:32 

22:28:33 

HH/VV 

25234 

130 

A 

27-Dec-06 

22:26:24 

22:26:39 

22:26:32 

22:25:33 

22:28:34 
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Table  2:  Statistics  for  validated  ship  signatures. 


Polarization 

= 

HH 

Total  #AIS  records  within  ASAR  coverage 

= 

2038 

#  AIS  records  with  IMO  <  0 

= 

0 

#  AIS  records  with  IMO  =  0 

= 

53 

#  AIS  records  with  IMO  =  1 

= 

0 

#  AIS  ships  landmasked 

= 

944 

#  AIS  ships  portmasked 

= 

174 

#  AIS  ships  at  image  edge 

= 

19 

#  AIS  ships  with  airballs 

= 

31 

#  AIS  ships  with  multitargs  detected  in  target  subscene 

= 

25 

#  AIS  ships  with  multiple  AIS  hits  at  same  target  location 

= 

5 

#  AIS  ships  with  no  significant  target  signature 

= 

78 

#  counts  where  no  AIS  to  ISR  match  occurred 

= 

20 

#  aistargtotalrcs<=0 

= 

0 

#  ISR  lengths  found  <=  0  meters 

= 

0 

#  AIS  ships  with  signiftargsig  >  search  radius  of  400m 

= 

17 

#  AIS  ships  with  signiftargsig  <=  search  radius  of  400m 

= 

672 

Polarization 

= 

w 

Total  #AIS  records  within  ASAR  coverage 

= 

1921 

#  AIS  records  with  IMO  <  0 

= 

0 

#  AIS  records  with  IMO  =  0 

= 

38 

#  AIS  records  with  IMO  =  1 

= 

0 

#  AIS  ships  landmasked 

= 

862 

#  AIS  ships  portmasked 

= 

164 

#  AIS  ships  at  image  edge 

= 

19 

#  AIS  ships  with  airballs 

= 

30 

#  AIS  ships  with  multitargs  detected  in  target  subscene 

= 

22 

#  AIS  ships  with  multiple  AIS  hits  at  same  target  location 

= 

3 

#  AIS  ships  with  no  significant  target  signature 

= 

84 

#  counts  where  no  AIS  to  ISR  match  occurred 

= 

30 

#  aistargtotalrcs<=0 

= 

0 

#  ISR  lengths  found  <=  0  meters 

= 

1 

#  AIS  ships  with  signiftargsig  >  search  radius  of  400m 

= 

10 

#  AIS  ships  with  signiftargsig  <=  search  radius  of 400m 

= 

658 

Polarization 

= 

HV 

Total  #AIS  records  within  ASAR  coverage 

= 

1083 

#  AIS  records  with  IMO  <  0 

= 

0 

#  AIS  records  with  IMO  =  0 

= 

28 

#  AIS  records  with  IMO  =  1 

= 

0 

#  AIS  ships  landmasked 

= 

506 

#  AIS  ships  portmasked 

= 

86 

#  AIS  ships  at  image  edge 

= 

6 

#  AIS  ships  with  airballs 

= 

19 

#  AIS  ships  with  multitargs  detected  in  target  subscene 

= 

13 

#  AIS  ships  with  multiple  AIS  hits  at  same  target  location 

= 

2 

#  AIS  ships  with  no  significant  target  signature 

= 

15 

#  counts  where  no  AIS  to  ISR  match  occurred 

= 

11 

#  aistargtotalrcs<=0 

= 

0 

#  ISR  lengths  found  <=  0  meters 

= 

0 

#  AIS  ships  with  signiftargsig  >  search  radius  of  400m 

= 

10 

#  AIS  ships  with  signiftargsig  <=  search  radius  of  400m 

= 

387 

Polarization 

= 

VH 

Total  #AIS  records  within  ASAR  coverage 

= 

966 

#  AIS  records  with  IMO  <  0 

= 

0 

#  AIS  records  with  IMO  =  0 

= 

13 

#  AIS  records  with  IMO  =  1 

= 

0 

#  AIS  ships  landmasked 

= 

424 

#  AIS  ships  portmasked 

= 

76 

#  AIS  ships  at  image  edge 

= 

6 

#  AIS  ships  with  airballs 

= 

18 

#  AIS  ships  with  multitargs  detected  in  target  subscene 

= 

6 

#  AIS  ships  with  multiple  AIS  hits  at  same  target  location 

= 

0 

#  AIS  ships  with  no  significant  target  signature 

= 

6 

#  counts  where  no  AIS  to  ISR  match  occurred 

= 

18 

#  aistargtotalrcs<=0 

= 

0 

#  ISR  lengths  found  <=  0  meters 

= 

0 

#  AIS  ships  with  signiftargsig  >  search  radius  of  400m 

= 

3 

#  AIS  ships  with  signiftargsig  <=  search  radius  of 400m 

= 

396 
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Dover, Gibraltar:  Error  Rate  =  10.4631% 


Figure  1:  Scatterplot  of  AISLive-reported  and  ISR  Database-validated  ship  lengths. 


Dover, Gibraltar:HH-Pol,  SearchRadius=400  (m),  n=672,  Mean=1 48.6044,  STD=77.1219 
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Dover, Gibraltar:VV-Pol,  SearchRadius=400  (m),  n=658,  Mean=1 52.3078,  STD=74.1941 
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Distance  Between  AlS-Predicted  Ship  Location  and  Maximum  o°  Location  of  Target  Signature  [m] 


Figure  2:  Histograms  of  the  distance  between  the  AlS-predicted  ship  position  and  the  SAR 
signature  position  for  each  linear  polarization. 
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ISR  Ship  Length  [m]  ISR  Ship  Length  [m] 


Dover, Gibraltar:  HV-Pol,  N=387p=1 45.3596, 0=64.7397 
25 1 - T -  - 1 - 1 - 1 - 1 - 
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Figure  3:  Histograms  of  ISR  Database  validated  ship  lengths  for  each  linear  polarization. 


DRDC  Ottawa  TM  2008-005 


7 


Dover, Gibraltar:  HH-Pol,  N  =672,n=42.97,  o=5.6587 
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Figure  4:  Histograms  of  total  RCS  for  the  validated  ships  for  each  linear  polarization. 
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Gibraltar,  2007-01-03,  IS1,  VV-pol  Gibraltar,  2007-01-03,  IS1,  VH-pol 


Figure  5:  An  example  ofVV  and  VH  images  of  ships  at  low  incidence  angles. 


ASA-APP-1  P, 2006-344-21 .25:02.328659,  Asc.Lstart  :3555,Lstop :  3612, MDSUVV  (Blue), MDS2=VH  (Green) 


15  16  17  18  19  20  21  22 

9jnc  [deg] 


ASA-APP-1  P, 2007-010-21  -.50:32.929413, Asc.Lstart :  2102,Lstop :  2169.MDSUHV  (Blue),MDS2=HH  (Green) 


Figure  6:  Range-ward  transects  of  normalized  radar  cross  section  through  ship  targets  at  near 

(left)  and  far  (right)  incidence  angles. 
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3  Ocean  clutter 


To  better  understand  the  clutter  against  which  the  ship  signature  is  measured,  the  clutter  estimates 
derived  during  the  ship  signature  analysis  were  plotted  as  a  function  of  incidence  angle  (0inc)  for 
co-polarizations  in  Figure  7  and  for  cross-polarizations  in  Figure  8.  Also  plotted  are  Envisat 
ASAR  noise  floor  estimates,  expressed  as  noise-equivalent-sigma-zero  [5]. 

These  plots  illustrate  the  following:  the  clutter  level  decreases  with  increasing  incidence  angle  for 
co-polarization,  sometimes  reaching  the  instrument  noise  floor;  VV  clutter  levels  are  generally 
higher  than  HH  clutter  levels;  the  cross-polarization  clutter  is  normally  at  the  noise  floor,  except 
for  the  smallest  incidence  angles  in  the  data  set  at  which  there  could  be  as  much  as  a  few  dB  of 
clutter  signal. 

In  ship  detection  systems,  the  background  clutter  level  is  generally  modelled  as  K-distributed 
noise.  We  represent  the  clutter  variability  by  plotting  in  Figure  9  the  mean-squared  to  variance 
ratio  (MSVR),  a  measure  of  the  number  of  statistically  independent  looks  for  a  Gaussian  scene,  as 
a  function  of  the  incidence  angle  and  polarization.  The  upper  limit  of  this  parameter  should  be 
close  to  the  number  of  independent  looks,  which  is  range-dependent  for  Envisat  ASAR  since  the 
ground  range  resolution  is  held  more-or-less  constant  with  range  and  the  extra  available  spatial 
resolution  is  converted  to  radiometric  resolution.  Also  plotted  are  nominal  values  of  the 
equivalent  number  of  looks  (ENL)  for  each  beam  mode  [5]. 

Smaller  values  of  MSVR  are  related  to  higher  clutter  variability,  indicating  non-Gaussian 
behaviour  that  may  be  characterized  through  use  of  a  K-distributed  clutter  model.  One  approach 
to  estimating  the  K-distribution  order  parameter  uses  the  observed  MSVR  and  the  ENL  [6]. 
Histograms  of  the  order  parameter  are  shown  in  Figure  10.  For  co-polarization,  66.5%  of  the 
cases  have  a  K-distribution  order  parameter  that  is  smaller  than  100,  representing  cases  with  a 
departure  from  Gaussian  statistics.  That  is,  the  co-polarization  clutter  is  generally  non-Gaussian 
in  nature.  On  the  other  hand,  for  cross-polarization,  only  18.6%  of  the  cases  have  a  K- 
distribution  order  parameter  that  is  smaller  than  100.  That  is,  the  cross-polarization  clutter  is 
generally  Gaussian  in  nature.  Indeed,  the  ASAR  cross-polarization  clutter  is  generally  just  the 
instrument  noise  floor. 
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Dover,  Gibraltar:  Envisat  ASAR,  HH-pol,  All  Ships,  N=672 


Figure  7:  Mean  co-polarization  ocean  clutter  as  a  function  of  incidence  angle  for  each  ship 

analyzed. 


Dover,  Gibraltar:  Envisat  ASAR,  HV-pol,  All  Ships,  N=387 
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Dover,  Gibraltar:  Envisat  ASAR,  VH-pol,  All  Ships,  N=396 


Figure  8:  Mean  cross-polarization  ocean  clutter  as  a  function  of  incidence  angle  for  each  ship 

analyzed. 
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Dover, Gibraltar:  HH-Pol,  SearchRadius=400  (m),  n=654  Dover, Gibraltar:  VV-Pol,  SearchRadius=400  (m),  n=633 


Figure  9:  Mean-squared  to  Variance  Ratio  as  a  function  of  incidence  angle  for  the  ocean  clutter 
associated  with  each  ship  analyzed  for  each  linear  polarization.  The  red  lines  indicate  the  ENL 

for  each  ASAR  beam  mode. 
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Dover, Gibraltar:  HH-Pol,  SearchRadius=400  (m),  n=619,  NK=397 
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Figure  10:  Histograms  of  the  K-distribution  order  parameter  for  the  ocean  clutter  associated 
with  each  ship  analyzed  and  for  each  linear  polarization. 
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4  Target  RCS  and  contrast  metrics 


In  developing  target  contrast  metrics,  we  considered  three  target  RCS  metrics:  Total  RCS,  Peak 
Region  RCS,  and  Segmented  Region  RCS.  All  are  based  upon  a  64  by  64  pixels  analysis  subscene 
that  is  centred  on  the  peak  G°  (i.e.,  normalized  radar  cross  section)  value.  The  16-bit  ESA 
products  considered  are  radiometrically  calibrated  and  provide  G°  directly.  Note  that  ship  target 
saturation  was  not  observed  in  any  of  the  datasets. 

Total  RCS  is  the  integration  of  <j°  over  the  entire  target  subscene  with  the  mean  clutter,  based 
upon  clutter  estimates  from  the  four  1 6  by  16  pixels  comers  of  the  analysis  subscene,  removed. 

Peak  Region  RCS  (PRrcs  )  is  the  integration  of  G°  over  the  peak  region  with  the  mean  clutter 

removed.  The  peak  region  is  defined  to  be  twice  the  3  dB  width  of  the  target  signature  centered 
on  the  maximum  G°  location  in  the  target  subscene  (this  is  a  standard  Envisat  ASAR  calibration 
team  definition). 

Segmented  Region  RCS  (  SRrcs  )  is  the  integration  of  G°  over  a  segmented  target  region  with  the 

mean  clutter  removed.  A  pixel  within  the  target  subscene  is  included  in  the  segmented  target 
region  if: 

cr°  >  mean  (g° dut )  +  5  std  (g° dut )  , 

where  mean  and  std  are  the  average  value  and  the  standard  deviation  operations.  The  segmented 
region  can  be  “cleaned-up”  by  using  sequential  erosion  and  dilation  morphological  operations, 
which  is  useful  for  target  length  estimation,  but  generally  has  little  impact  on  the  estimated 


Some  possible  contrast  metrics  include:  Peak-to-Clutter  Ratio,  Peak  Region  RCS-to-Clutter 
Ratio,  and  the  Segmented  Region  RCS-to-Clutter  Ratio. 

Peak-to-Clutter  Ratio  is  defined  as  G°  /  mean  (g°  , , ) ,  where  G°  is  the  target  subscene 
maximum  <j°  value. 

Peak  Region  RCS-to-Clutter  Ratio  is  defined  as  [  PRrcs  /  mean  ( G°dui )] /  PRarea  where  PRarea  is 
the  Peak  Region  area  in  square  meters; 

Segmented  Region  RCS-to-Clutter  Ratio  is  equal  to  [  SR  RCS  /  mean  (g° du, )]/  SRarea  where  SRarea 
is  the  Segmented  Region  area  in  square  meters. 

Scatterplots  of  the  peak  region  RCS  and  the  total  RCS  are  shown  in  Figure  11.  In  general,  the 
peak  region  RCS  underestimates  the  total  RCS.  The  peak  region  RCS  is  not  appropriate  if  the 
signature  contains  several  bright  components  since  only  the  largest  is  measured.  Discrepancies  of 
up  to  15  dB  are  noted. 
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Scatterplots  of  the  segmented  region  RCS  and  the  total  RCS  are  shown  in  Figure  12.  These  RCS 
measures  are  generally  in  good  agreement,  especially  for  larger  RCS  values.  For  contrast 
measures  we  will  generally  focus  on  using  the  segmented  region  RCS  since  a  target  area  estimate 
is  required  for  the  contrast  metric. 

Scatterplots  of  the  peak-to-clutter  ratio  and  the  peak  region  RCS-to-clutter  ratio  are  shown  in 
Figure  13.  These  contrast  metrics  are  well-correlated  for  larger  contrast  values. 

Scatterplots  of  the  peak-to-clutter  ratio  and  the  segmented  region  RCS-to-clutter  ratio  are  shown 
in  Figure  14.  These  contrast  metrics  are  reasonably  well  correlated,  but  we  feel  that  the 
segmented  region  RCS-to-clutter  ratio  is  a  better  measure  since  it  includes  an  RCS  that  is  closer 
to  the  total  RCS.  For  completeness,  all  of  the  contrast  metrics  are  included  in  subsequent 
analysis. 


Dover,  Gibraltar:  Envisat  ASAR,  HH-pol,  N=672  Dover,  Gibraltar:  Envisat  ASAR,  VV-pol,  N=658 


Dover,  Gibraltar:  Envisat  ASAR,  HV-pol,  N=387  Dover,  Gibraltar:  Envisat  ASAR,  VH-pol,  N=396 


Figure  11:  Scatterplots  of peak  region  RCS  and  total  RCS  for  each  linear  polarization. 
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Figure  12:  Scatterplots  of  segmented  region  RCS  and  total  RCS  for  each  linear  polarization. 


Dover,  Gibraltar:  Envisat  ASAR,  HH-pol,  N=672 


Dover,  Gibraltar:  Envisat  ASAR,  HV-pol,  N=387 


Dover,  Gibraltar:  Envisat  ASAR,  VV-pol,  N=658 


Dover,  Gibraltar:  Envisat  ASAR,  VH-pol,  N=396 


Figure  13:  Scatterplots  of peak-to-clutter  ratio  and  peak  region  RCS-to-clutter  ratio  for  each 

linear  polarization. 
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Dover,  Gibraltar:  Envisat  ASAR,  HH-pol,  N=672 


HH-pol  Segmented  Region  RCS-to-Clut  Ratio  [dB] 
Dover,  Gibraltar:  Envisat  ASAR,  HV-pol,  N=387 


HV-pol  Segmented  Region  RCS-to-Clut  Ratio  [dB] 


Dover,  Gibraltar:  Envisat  ASAR,  HH-pol,  N=658 


VV-pol  Segmented  Region  RCS-to-Clut  Ratio  [dB] 
Dover,  Gibraltar:  Envisat  ASAR,  VH-pol,  N=396 


VH-pol  Segmented  Region  RCS-to-Clut  Ratio  [dB] 


Figure  14:  Scatterplots  of peak-to-clutter  ratio  and  segmented  region  RCS-to-clutter  ratio  for 

each  linear  polarization. 


4.1  Total  RCS 

Scatterplots  of  the  total  RCS  versus  the  ISR  Database  validated  ship  length  for  each  polarization 
are  shown  in  Figure  15.  In  each  case,  a  simple  regression  of  the  RCS  as  a  function  of  the  ship 
length  has  been  carried  out  to  show  the  general  trend  in  the  data;  the  model  coefficients  are 
included  in  the  title  of  each  sub-plot.  For  comparison,  previous  such  model  fits  to  the  R-l  Fine 
mode  [2]  and  SCNB  [3]  data  are  included  on  each  plot.  We  see  that  the  co-polarization  data  are 
very  similar  to  the  R-l  data  acquired  at  HH  polarization,  while  the  cross-polarization  data  are 
around  10  dB  smaller,  as  observed  previously  with  Environment  Canada  (EC)  CV-580  C-band 
polarimetric  SAR  imagery  of  known  ships  [7]. 

Scatterplots  of  total  RCS  for  co-  and  cross-polarization  data  (i.e.,  combined  HH  and  VV  data,  and 
combined  HV  and  VH  data)  are  shown  in  Figure  16,  along  with  the  combined  model  fits.  At  this 
stage,  we  don’t  observe  any  obvious  differences  between  the  co-polarization  channels  (i.e., 
between  HH  and  VV)  or  between  the  cross-polarization  channels  (i.e.,  between  VH  and  HV). 

Table  3  summarizes  the  ASAR  total  RCS  regression  fit  parameters  where  a  is  a  function  of  the 
ship  length  (L).  The  R-l  total  RCS  regression  fit  parameters  are  also  included  for  both  the  Fine 
and  SCNB  mode  datasets. 
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Figure  17  shows  the  total  RCS  regression  fits  for  the  R-l  and  ASAR  data  sets  on  the  left,  and  the 
total  RCS  regression  fits  normalized  by  the  RADARS  AT  Fine  mode  fit  on  the  right.  These  plots 
illustrate  that  the  co-polarization  results  agree  among  the  various  regression  fits  to  within  2  dB, 
for  all  ship  lengths.  On  the  other  hand,  the  cross-polarization  ship  RCS  is  about  10  dB  smaller 
than  the  co-polarization  ship  RCS. 


Dover, Gibraltar:  HH-Pol,  SearchRadius=400  (m),  N=672,o=2.0627L**1.8649 


Dover, Gibraltar:  VV-Pol,  SearchRadius=400  (m),  N=658,o=0.77938L**2.1183 


Dover, Gibraltar, HV-Pol,  SearchRadius=400  (m),  N=387,o=0.31276L"1.8084 


Dover, Gibraltar:  VH-Pol,  SearchRadius=400  (m),  N=396,u=0.14625L**1 .9446 


Figure  15:  Scatterplots  of  total  RCS  and  1SR  Database  validated  ship  length  for  each  linear 

polarization. 


18 


DRDC  Ottawa  TM  2008-005 


HCS|rf] 


Figure  16:  Scatterplots  of  total  RCS  and  ISR  Database  validated  ship  length  for  co-polarization 

(left)  and  cross-polarization  (right)  data. 


Table  3:  Total  RCS  regression  fit  parameters. 


Regression  Equation 

#  Ships 

ASAR 

oHH  =  2.0627  *LL8649 

672 

66 

ovv  =  0.77938  *  L2'1183 

658 

66 

cv  =0.3 1276*  L1  8084 

387 

66 

oVH  =0.1 4625*  L1  9446 

396 

66 

<j  =  1.2432  *L‘ 9956 

Co— pot 

1330 

66 

(Jr  ,=  0.22009  *L‘  8707 

Cross— pot 

783 

R-l 

oF.  =  1.2499  *L‘  9798 

399 

66 

°scm  =  0.32425  *L2'2991 

1805 

DRDC  Ottawa  TM  2008-005 


19 


Model  Fits  vs  Ship  Length 


Model  Fit  Comparis 


Ship  Length 


Figure  17:  Total  RCS  regression  fits  for  the  R-l  and  ASAR  data  sets  (left)  and  the  total  RCS 
regression  fits  normalized  by  the  R-l  Fine  mode  fit  (right). 

4.2  RCS  and  contrast  metrics  versus  ship  length 

We  now  consider  the  extracted  ship  contrast  metrics  as  a  function  of  the  validated  ISR  ship 
length.  The  three  RCS  metrics  for  co-polarization  have  been  plotted  in  Figure  18  and  for  cross¬ 
polarization  in  Figure  19.  In  each  case,  the  R-l  Fine  and  SCNB  models  were  based  on  the  total 
RCS. 

The  three  contrast  metrics  for  co-polarization  are  shown  in  Figure  20  and  for  cross-polarization  in 
Figure  21.  Included  in  the  plots  are  linear  fits  to  the  observed  data  to  show  the  general  trend.  The 
dependence  on  ship  length  is  much  clearer  for  the  cross-polarization  data,  which  is  essentially 
measured  against  the  instrument  noise  floor.  It  is  apparent  that  the  cross-polarization  data  have 
much  less  variability. 
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Dover,  Gibraltar:  HH  &  VV  Pol,  N=1330,g=1  ,2432L**1 ,9990over,  Gibraltar:  HH  &  VV  Pol,  N=1330,g=1.1526L**1.874 


ISR  Database  Ship  Length  [m]  ISR  Database  Ship  Length  [m] 


Dover,  Gibraltar:  HH  &  VV  Pol,  N=1330,o=0.63878L**2.102 


ISR  Database  Ship  Length  [m] 


Figure  18:  Scatterplots  ofRCS  metrics  for  co-polarization  versus  validated  ship  length. 


Dover,  Gibraltar:  HV  &  VH  Pol,  N=783,o=0.22009L**1 ,87ID3bver,  Gibraltar:  HV  &  VH  Pol,  N=783,o=0.58861L**1.5265 


ISR  Database  Ship  Length  [m]  ISR  Database  Ship  Length  [m] 


Dover,  Gibraltar:  HV  &  VH  Pol,  N=783,o=0.13151L**1.9572 


ISR  Database  Ship  Length  [m] 


HV  &  VH  RCS  Measure 
HV  &  VH  RCS/ASAR  Fit 
HH  RCS/RSAT  Fine  Mode  Fit 
HH  RCS/RSAT  SCNB  Fit 


Figure  19:  Scatterplots  ofRCS  metrics  for  cross-polarization  versus  validated  ship  length. 
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Dover,  Gibraltar:  HH  &  VV  Pol,  N=1330,a=0.027551L**1D68©r,  Gibraltar:  HH  &  VV  Pol,  N=1 330, a=0.0052033L**1 .8228 


Figure  20:  Scatterplots  of  contrast  metrics  for  co-polarization  versus  validated  ship  length. 


Dover,  Gibraltar:  HV  &  VH  Pol,  N=783,a=0.1426L**1 .46EGbver<cGibraltar:  HV  &  VH  Pol,  N=783,o=0.031866L**1 .5104 


^  ISR  Database  Ship  Length  [m] 


ISR  Database  Ship  Length  [m] 


Dover, Gibraltar:  HV  &  VH  Pol,  N=783,a=0.44503L**0.83262 


HV  &  VH  Contrast  Measure 
Fit 


Figure  21:  Scatterplots  of  contrast  metrics  for  cross-polarization  versus  validated  ship  length. 
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4.3  RCS  and  contrast  metrics  versus  incidence  angle 


We  now  consider  the  extracted  RCS  and  contrast  metrics  as  a  function  of  the  local  incidence 
angle.  The  metrics  for  each  polarization  are  shown  in  Figure  22  through  Figure  25.  The  data  gap 
around  37°  incidence  angle  is  due  to  1S5  mode  not  being  available  at  the  time  of  the  data 
acquisition  campaign.  Again,  linear  fits  to  the  data  are  included  in  each  case  so  that  the  general 
trend  is  apparent. 

Considering  the  segmented  region  contrast  (the  lower-right  plot  in  each  case),  the  target  contrast 
is  seen  to  improve  with  increasing  incidence  angle  for  co-polarization.  This  reflects  the  fact  that 
the  co-polarization  clutter  decreases  with  incidence  angle,  which  drives  recommendations  to  use 
larger  incidence  angles  for  ship  detection.  On  the  other  hand,  the  target  contrast  is  essentially 
independent  of  incidence  angle  for  cross-polarization.  This  is  due  to  the  clutter  being  noise  floor 
limited  for  cross-polarization.  Recall  that  the  noise  floor  is  essentially  independent  of  incidence 
angle. 

The  data  fdtering  algorithm  rejects  low  contrast  (i.e.  ,  <  5  dB  peak-to-clutter  ratio)  targets,  which 
is  evident  as  a  lower  limit  in  the  contrast  measures.  Therefore,  the  best  fit  straight  line  that  is 
plotted  to  show  the  general  dependence  on  incidence  angle  has  a  slope  that  is  biased  low. 

The  segmented  region  RCS-to-clutter  ratio  as  a  function  of  incidence  angle  is  shown  in  Figure  26 
for  both  co-polarization  (on  the  left)  and  cross-polarization  (on  the  right).  The  data  have  been 
broken  down  by  ship  length  population  with  the  total  population  at  the  top,  then  the  1/3  smallest 
ships  (L  <  100  m),  the  middle  1/3  ships  (100  m  <  L  <  160  m),  and  the  1/3  largest  ships  (L  >160 
m).  The  straight  line  fits,  which  are  summarized  in  Figure  27,  show  how  the  observed  contrast 
depends  on  the  incidence  angle,  ship  length,  and  polarization  (for  co-polarization  and  cross¬ 
polarization,  at  least).  These  plots  illustrate  how  ships  are  more  detectable  for  co-polarization 
data  at  larger  incidence  angles,  and  are  more  detectable  for  cross-polarization  data  at  smaller 
incidence  angles.  The  break-point  between  cross-polarization  and  co-polarization  providing  the 
best  contrast  depends  on  the  ship  length,  but  is  seen  to  be  around  33°  (i.e.,  it  is  best  to  use  cross¬ 
polarization  for  incidence  angles  smaller  than  33°).  This  dependence  on  incidence  angle  reflects 
the  decreasing  clutter  with  increasing  incidence  angle  for  co-polarization,  the  noise  floor  limited 
clutter  for  cross-polarization,  and  the  fact  that  the  cross-polarization  ship  target  RCS  is  about  10 
dB  smaller  than  the  co-polarization  ship  target  RCS. 
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Dover,  Gibraltar:  Envisat  ASAR,  HH-pol,  All  Ships,  N=673Dover,  Gibraltar:  Envisat  ASAR,  HH-pol,  All  Ships,  N=672 
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Figure  22:  RCS  and  contrast  metrics  as  a  function  of  incidence  angle  for  HH  polarization. 


Dover,  Gibraltar:  Envisat  ASAR,  VV-pol,  All  Ships,  N=65£Dover,  Gibraltar:  Envisat  ASAR,  VV-pol,  All  Ships,  N=658 
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Figure  23:  RCS  and  contrast  metrics  as  a  function  of  incidence  angle  for  VV  polarization. 
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Figure  24:  RCS  and  contrast  metrics  as  a  function  of  incidence  angle  for  HV  polarization. 
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Figure  25:  RCS  and  contrast  metrics  as  a  function  of  incidence  angle  for  VH  polarization. 
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Figure  26:  Segmented  region  RCS-to-clutter  ratio  for  co-polarization  (left)  and  cross- 
polarization  (right)  as  a  function  of  incidence  angle  for  (top-to-bottom)  all  ships,  smallest  ships, 
medium  ships,  and  longest  ships.  (Continued  on  next  page.) 
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Dover,  Gibraltar:  Envisat  ASAR,  HH  &  W  pols,  SL>1 60m,  N=525 


Dover,  Gibraltar:  Envisat  ASAR,  HV  8  VH  pols,  SL>160m,  N=291 


0|nc  (deg] 


Figure  26:  Concluded. 


Dover,  Gibraltar:  Envisat  ASAR,  Fit  Comparison 


Figure  27:  Summary >  of  contrast  as  a  function  of  incidence  angle  and  ship  length  for  both  co- 

polarization  and  cross-polarization. 


DRDC  Ottawa  TM  2008-005 


27 


4.4  AP  mode  comparisons 


Figure  28  shows  scatterplots  of  total  RCS  measurements  for  each  ship  in  the  available  AP  modes. 
There  are  three  such  plots  since  there  are  three  different  AP  modes  available  (i.e.,  HH/VV, 
HH/HV,  and  VV/VH). 

From  the  plots,  we  see  that  the  FIH/W  total  RCS  values  are  well  correlated.  Flowever,  the  co- 
polarization/cross-polarization  Total  RCS  values  do  not  appear  to  be  as  well  correlated.  Once 
again,  the  cross-polarization  total  RCS  is  about  10  dB  smaller  than  the  co-polarization  total  RCS. 


Dover,  Gibraltar:  Envisat  ASAR,  HH  vs  VV  pol,  N=285 


Dover,  Gibraltar:  Envisat  ASAR,  VV  vs  VH  pol,  N=352  Dover,  Gibraltar:  Envisat  ASAR,  HH  vs  HV  pol,  N=352 


Figure  28:  Scatterplots  of  total  RCS  for  each  AP  mode. 
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5  PDF  fits  to  the  normalized  total  RCS 


For  each  linear  polarization,  we  tested  various  Probability  Density  Function  (PDF)  models  to 
determine  if  they  can  describe  the  observed  variability  in  the  normalized  total  RCS  distribution 
including  Lognormal,  Gamma,  Weibull,  Rayleigh,  Ricean  and  Chi-Squared  [1].  The 
normalization  was  carried  out  by  dividing  the  observed  total  RCS  value  by  the  appropriate  fitted 
total  RCS  of  Table  3.  The  Kolomogorov-Smimov  (K-S)  goodness-of-fit  test  was  used  to 
determine  if  the  proposed  PDF  describes  the  observed  variability  (i.e.,  passes  or  fails). 

From  Table  3,  we  see  that  the  Lognormal  PDF  describes  the  normalized  RCS  variability  for  both 
HH  and  VV  polarizations,  whereas  the  Gamma  and  Weibull  PDFs  are  appropriate  for  HV 
polarization,  and  Lognormal  and  Gamma  PDFs  are  appropriate  for  VF1  polarization.  Figure  29 
provides  both  the  PDF  and  the  corresponding  Empirical  Cumulative  Density  Function  (ECDF) 
that  passed  the  K-S  test  for  each  linear  polarization.  Figure  30  shows  both  the  PDF  and  ECDF 
fits  that  passed  the  K-S  test  for  the  combined  co-polarization  and  cross-polarization  data  sets. 
Table  4  summarizes  the  results  and  provides  the  model  fit  parameters1  along  with  the  10%  and 
90%  cumulative  probability  values  for  the  fitted  PDFs.  The  results  indicate  that  the  co¬ 
polarization  data  are  less  dispersed  than  the  cross-polarization  data.  Furthermore,  the  provided 
PDFs  could  be  used  as  inputs  to  ship  detection  software  to  guide  the  probability  of  missed 
detection,  or  could  be  used  to  improve  ship  detectability  prediction  for  future  SAR  missions. 

Table  4:  Model  fit  results  for  each  linear  polarization  and  combined  polarizations. 


Pol 

Model  Fit 

A 

B 

10% 

90% 

HH 

Lognormal 

5.7560e-15 

1.0586 

0.26 

3.88 

VV 

Lognormal 

-5.3544e-15 

1.0758 

0.25 

3.97 

HV 

Gamma 

1.7189 

0.79979 

0.32 

2.78 

HV 

Weibull 

1.4999 

1.3352 

0.39 

4.16 

VH 

Lognormal 

2.7560e-15 

0.70430 

0.41 

2.47 

VH 

Gamma 

2.1870 

0.58462 

0.37 

2.44 

Co-Pol 

Lognormal2 

-6.1020e-3 

1.0661 

0.25 

3.89 

Cross-Pol 

Gamma 

1.9049 

0.69799 

0.34 

2.62 

1  In  Table  4  and  the  subsequent  plots,  the  A  and  B  parameters  correspond  respectively  to  the  Shape  and 
Scale  parameters  for  the  Gamma  distribution,  the  Scale  and  Shape  parameters  for  the  Weibull  distribution, 
and  mu  and  sigma  for  the  Lognormal  distribution.  The  PDFs  are  as  defined  in  Matlab " . 

2  The  two  largest  normalized  total  RCS  values  were  filtered  from  the  dataset  in  order  for  the  model  to  pass 
the  K-S  test. 
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Density 


Empirical  CDFs,  K-S  Hypothesis  Accepted,  N=672  P=0.13123,  KSSTAT=0.044803 


AIS/ASAR  Dover, Gibraltar,  HH-pol,  All  Shiplengths,  N=672,p=1.8767, 0=3.8613,  (i2/o2=0.23622,  =9.0142,  y2=112.6541 


Wrn  i 


—  Empirical  PDF 

—  Bar  Graph  Empirical  PDF 

—  logn  Fit,  Mu=5.756e-1 5,  Sigma=1 .0586 


95%  Confidence  Interval:  LB  =  -0.080182, 1 .0049 
95%  Confidence  interval:  U  B  =  0.0801 82,:  1 . 1 1 84 
K-S  Hypothesis  Accepted,  P-0.1 31 23,  KSSTAT=0.044803_ 

1 0%  Normalized  TotaIRCS  postion  on  PDF  =  0.26 
20%  Normalized  TotaIRCS  postion  on  PDF  =  0.41 
80%  Normalized  TotaIRCS  postion  on  PDF  =  2.44 
90%  Normalized  TotaIRCS  postion  on  PDF  =  3!88  :  _ 

99%  Normalized  TotaIRCS  postion  on  PDF  =  1 1 .71 


T  otal  RCS/Model  Fit 


AIS/ASAR  Dover, Gibraltar,  VV-pol,  All  Shiplengths,  N=658,p=1 .9343,  o=4.0812,  p2/o2= 0.22465,  ^=7.9188,  y2=83.2349 


— —  Empirical  PDF 

- Bar  Graph  Empirical  PDF 

- logn  Fit,  Mu=-5.3544e-1 5,  Sigma=1 .0758 


95%  Confidence  Interval  :LB  =  -0.082351 , 1:.0207 
95%  Confidence  Interval  :UB  =  0.082351 , 1.1373 
.  K-S.  Hypothesis  Accepted,  P=0.095975,  KSSTAT=.Q.b.47805 

1 0%  Normalized  TotaIRCS  postion  on  PDF  :=  0.25 
20%  Normalized  TotaIRCS  postion  on  PDF  :=  0.4 
80%  Normalized  TotaIRCS  postion  on  PDF  :=  2.47 
90%  Normalized  TotaIRCS  postion  on  PDF  :=  3.97  "I 
99%  Normalized  TotaIRCS  postion  on  PDF  :=  12.16  : 


20  25  30  35  40  45  50 

TotaIRCS/ModelFit 


Empirical  CDFs,  K-S  Hypothesis  Accepted,  N=658  P=0.095975,  KSSTAT=0.047805 


AIS/ASAR  Dover, Gibraltar,  HV-pol,  All  Shiplengths,  N=387,p=1.3747,  a=1. 1 291,  (i2/o2=1. 4825,  yfl  =2.9966,  y2=1 8.5607  Empirical  CDFs,  K-S  Hypothesis  Accepted,  N=387  P=0.30534,  KSSTAT=0.048918 


Figure  29:  PDF  (left)  and  ECDF  (right)  for  each  linear  polarization  that  passed  the  K-S 
goodness-of-fit  test.  (Continued  on  next  page.) 
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TotaIRCS/ModelFit 


Empirical  CDFs,  K-S  Hypothesis  Accepted,  N=387  P=0.08788,  KSSTAT=0.063136 


AIS/ASAR  Dover, Gibraltar,  VH-pol,  All  Shiplengths,  N=396,p=1 .2785,  c=1 .1209,  p2/o2=1 .301 1 ,  y,  =5.6773,  y2=63.4447 


TotaIRCS/ModelFit 


Empirical  CDFs,  K-S  Hypothesis  Accepted,  N=396  P=0.43768,  KSSTAT=0.043369 


AIS/ASAR  Dover, Gibraltar,  VH-pol,  All  Shiplengths,  N=396,p=1.2785,  a=1.1209,  p2/cr2=1 .301 1 ,  yt =5.6773,  y2=63.4447 


TotaIRCS/ModelFit 


Empirical  CDFs,  K-S  Hypothesis  Accepted,  N=396  P=0.092053,  KSSTAT=0.061954 


Figure  29:  Concluded. 
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AIS/ASAR  Dover, Gibraltar,  HHVV-pol,  All  Shiplengths,  N=1 330,  p=1 .9229,  a=4.0017,  |i2/o2=0.23091 ,  y  =8.2282,  y2=91. 971 8 


T  otaIRCS/Model  Fit 


AIS/ASAR  Dover, Gibraltar,  HHVV-pol,  All  Shiplengths,  N=1 328,  (1=1.8392, 0=3.3716,  p2/o2=0.29756, 7,  =7.0696, 


- Empirical  PDF 

- Bar  Graph  Empirical  PDF 

- logn  Fit,  Mu=-0.0061 02,  Sigma=1.0661 


95%  Confidence  Interval  LB  =  -0-063491 , 1.027 
95%  Confidence  Interval  UB  =  0.051 287, 1 .1 082 
K-S  Hypothesis  Accepted,  P=0.Q60668,  KSSTAT=0.0361 56 

1 0%  Normalized  TotaIRCS  postion  on  PDF  :=  0.25 
20%  Normalized  TotaIRCS  postion  on  PDF  :=  0.41 
80%  Normalized  TotaIRCS  postion  on  PDF  :=  2.44 
90%  Normalized  TotaIRCS  postion  on  PDF :='  3.89 
99%  Normalized  TotaIRCS  postion  on  PDF  :=  1 1 .8 


20  25  30  35  40  45 

TotaIRCS/ModelFit 


y2=71.1816 


AIS/ASAR  Dover, Gibraltar,  HVVH-pol,  All  Shiplengths,  N=783,  p=1.3296,  o=1.1379,  ji2/ct2=1  .3654,  y,  =4.2261,  y2=37.4323 


TotaIRCS/ModelFit 


Figure  30:  Combined  co-polarization  and  cross-polarization  PDF  (left)  and  ECDF  (right).  The 
original  co-polarization  data  set  (Is'  row)  did  not  pass  the  K-S  test;  however,  it  does  pass  the  K-S 
test  (2nd  row)  after  removing  the  two  largest  normalized  total  RCS  samples  from  the  data  set. 
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6  Target  signature  length  estimates 


Although  there  are  many  approaches  to  ship  signature  length  estimation,  in  this  work  we  have 
focussed  on  the  results  from  VUSAR  [8],  a  tool  that  is  in  routine  use  at  DRDC  Ottawa  and  has 
been  extended  for  ship  signature  analysis.  VUSAR  segments  the  ship  signature  from  the 
background  clutter  by  using  a  thresholding  operation,  followed  by  a  sequence  of  morphological 
operations  to  better  define  the  signature’s  shape.  Derived  ship  length  estimates  are  expected  to  be 
representative  of  results  from  other  ship  signature  analysis  tools.  Figure  3 1  shows  scatterplots  of 
ISR  database  validated  ship  length  and  the  VUSAR- measured  ship  length  estimates  for  each 
linear  polarization.  We  note  that  there  is  fairly  good  correlation  for  all  polarizations,  but  that 
there  are  occurrences  of  ship  length  under  and  over  estimation.  The  occurrences  of  these 
conditions  will  be  considered  in  greater  detail. 

Figure  32  provides  histograms  of  the  error  measurement  s  =  AL/L  where  AL  =  (VUSAR 
measured  ship  signature  length  -  ISR  database  validated  ship  length)  and  L  is  the  ISR  database 
validated  ship  length.  The  red  vertical  lines  to  the  left  and  right  of  zero  in  each  plot  represent  the 
30%  error  tolerance  locations  that  were  considered  previously.  The  rates  of  occurrence  are 
summarized  in  Table  5.  The  “e  <  -0.3”  column  represents  the  percentage  of  ship  length  estimates 
below  the  30%  error  tolerance,  while  the  “s  >  0.3”  represents  the  percentage  of  ship  length 
estimates  above  the  30%  error  tolerance.  For  co-polarization  we  see  a  considerably  larger  rate  of 
underestimation  of  ship  length  (approximately  27%)  than  for  cross-polarization  (approximately 
6%).  Presumably  this  arises  since  the  clutter  levels  are  higher,  often  masking  part  of  the  ship 
signature.  On  the  other  hand,  for  cross-polarization  we  see  a  considerably  larger  rate  of 
overestimation  of  ship  length  (approximately  22%)  than  for  co-polarization  (approximately  12%). 
This  result  is  somewhat  surprising  since  we  expect  that  the  degree  of  azimuth  signature  smearing 
due  to  ship  motion  should  be  independent  of  polarization.  However,  it  is  possible  that  the 
generally  high  contrast  ratio  for  cross-polarization  might  more  clearly  show  signature  sidelobes  or 
other  imaging  artefacts  that  are  being  measured  as  part  of  the  ship  length. 

Table  5:  Breakdown  of  all  ship  length  estimates. 


All  Ships 

Pol 

N 

£  <  -0.3 

[%1 

£>0.3 

[%] 

HH 

654 

26.8 

13.3 

VV 

633 

28.9 

10.9 

HV 

373 

6.7 

23.6 

VH 

392 

6.1 

20.7 

In  the  balance  of  this  section,  we  investigate  dependencies  of  target  signature  length  estimates  on 
incidence  angle,  target-to-clutter  ratio,  and  ship  aspect  angle. 

Figure  33  shows  scatterplots  of  ISR  database  validated  ship  length  and  the  VUSAR  measured 
ship  length  for  each  linear  polarization.  The  red  points  represent  ships  with  9inc  <  25°  and  the 
blue  points  represent  ships  with  9inc  >  25°.  For  co-polarization,  there  are  more  underestimates  of 
ship  length  at  lower  incidence  angles,  as  observed  by  the  large  number  of  red  points  that  appear 
below  the  1:1  fit  line.  For  cross-polarization,  there  are  very  few  underestimates  observed  at  all 
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incidence  angles.  There  are  more  overestimates  of  ship  length  for  cross-polarization.  Figure  34 
provides  the  AL/L  histogram  for  ships  with  0inc  <  25°,  while  Figure  35  provides  the  AL/L 
histogram  for  ships  with  9inc  >  25°.  Table  6  summarizes  the  results. 


Table  6:  Breakdown  of  ship  length  estimates  by  incidence  angle. 


All  Ships 

0i„c  <  25° 

0i„c  >  25° 

Pol 

N 

£  <  -0.3 

[%1 

£>0.3 

[%] 

N 

£  <-  0.3 

[%1 

£>0.3 

[%] 

N 

£  <  -0.3 

[%1 

£>0.3 

|%] 

HH 

654 

26.8 

13.3 

171 

64.9 

5.8 

483 

13.2 

17.0 

VV 

633 

28.9 

10.9 

174 

64.4 

4.0 

459 

14.8 

13.7 

HV 

373 

6.7 

23.6 

117 

9.4 

30.8 

256 

5.9 

20.3 

VH 

392 

6.1 

20.7 

154 

7.1 

23.4 

238 

5.4 

23.3 

Figure  36  shows  scatterplots  of  ISR  database  validated  ship  length  and  the  VUSAR  measured 
ship  length  for  each  linear  polarization.  In  this  case,  the  red  points  represent  ships  with  TCR3  < 
10  dB  and  the  blue  points  represent  ships  with  TCR  >10  dB.  For  co-polarization,  a  significant 
number  of  ships  are  underestimated  at  low  TCRs,  as  observed  by  the  large  number  of  red  samples 
that  fall  below  the  1:1  fit  line.  For  cross-polarization,  there  are  very  few  underestimates 
observed.  There  appears  to  be  more  overestimates  of  ship  length  for  cross-polarization.  Figure 
37  provides  the  AL/L  histogram  for  ships  with  TCR  <  10  dB,  while  Figure  38  provides  the  AL/L 
histogram  for  ships  with  TCR  >10  dB.  Table  7  summarizes  the  results. 


Table  7:  Breakdown  of  ship  length  estimates  by  TCR. 


All  Ships 

TCR  <  10  dB 

TCR  >  10  dB 

Pol 

N 

£  <  -0.3 

[%] 

£>0.3 

[%] 

N 

£  <  -0.3 

[%1 

£>0.3 

[%] 

N 

£  <  -0.3 

[%1 

£>0.3 

[%] 

HH 

654 

26.8 

13.3 

175 

63.4 

9.1 

479 

13.8 

16.3 

VV 

633 

28.9 

10.9 

157 

66.9 

9.6 

476 

16.2 

10.5 

HV 

373 

6.7 

23.6 

29 

17.2 

41.4 

344 

5.8 

22.1 

VH 

392 

6.1 

20.7 

21 

28.6 

33.3 

371 

4.9 

20.2 

Figure  39  shows  scatterplots  of  ISR  database  validated  ship  length  and  the  VUSAR  measured 
ship  length  for  each  linear  polarization.  In  this  case,  the  red  points  represent  azimuth  travelling 
ships  (i.e.,  within  ±45°  of  travelling  along  the  azimuth  direction,  parallel  or  anti-parallel  to  the 
satellite  track)  and  the  blue  points  represent  range  travelling  ships  (i.e.,  within  ±45°  of  travelling 
along  the  range  direction,  towards  or  away  from  the  satellite).  For  both  co-polarization  and  cross¬ 
polarization,  most  of  the  larger  overestimates  are  azimuth  travelling  ships.  Figure  40  provides  the 
AL/L  histogram  for  range  travelling  ships  and  Figure  4 1  provides  the  AL/L  histogram  for  azimuth 
travelling  ships.  Table  8  summarizes  the  results. 


3  Represented  by  the  segmented  region  RCS-to-clutter  ratio. 
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Table  8:  Breakdown  of  ship  length  estimates  by  ship  aspect  angle. 


All  Ships 

Range  Travelling 

Azimuth  Travelling 

Pol 

N 

£  <  -0.3 

[%1 

£>0.3 

[%] 

N 

£  <  -0.3 

[%1 

£>0.3 

[%] 

N 

£  <  -0.3 

[%1 

£>0.3 

[%] 

HH 

654 

26.8 

13.3 

444 

22.8 

9.9 

210 

33.3 

20.5 

VV 

633 

28.9 

10.9 

487 

27.7 

OO 

bo 

146 

32.2 

15.1 

HV 

373 

6.7 

23.6 

234 

7.7 

10.3 

139 

4.3 

48.9 

VH 

392 

6.1 

20.7 

318 

6.6 

14.8 

74 

4.1 

51.4 

Figure  31:  Scatterplots  ofVUSAR  ship  length  and  ISR  Database  validated  ship  length  for  each 

linear  polarization. 
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Number  of  Ships  Number  of  Ships 


Figure  32:  Histograms  of  AL/L  showing  the  30%  error  positions  (red  lines)  for  each  linear 

polarization. 
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Dover, Gibraltar:  HH-Pol,  SearchRadius=400  (m),  n=654 


Validated  ship  length  [m] 


Dover, Gibraltar:  VV-Pol,  SearchRadius=400  (m),  n=633 


Validated  ship  length  [m] 


Dover, Gibraltar:  HV-Pol,  SearchRadius=400  (m),  n=373 


Dover, Gibraltar:  VH-Pol,  SearchRadius=400  (m),  n=392 


Figure  33:  Scatterplots  ofVUSAR  ship  length  and  ISR  Database  validated  ship  length  for  each 
linear  polarization.  The  blue  points  indicate  ships  at  incidence  angles  >  25°  while  the  red  points 

indicate  ships  at  incidence  angles  <  25°. 
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Dover, Gibraltar:  HH-Pol,|i=-0.38266,  o=0.42092,  n=171,  %Under30=64.91 23,  %Over30=5.848 


Dover, Gibraltar:  VV-Pol^=-0.40512,  0=0.36264,  n=174,  %Under30=64.3678,  %0ver30=4.023 


Dover, Gibraltar:  VH-Pol,|i=0.20528,  a=0.58481,  n=154,  %Under30=7.1429,  %Over30=23.3766 


AL/L 


Figure  34:  Histograms  of  AL/L  for  ships  at  incidence  angles  <  25°  showing  the  30%  error 
positions  (red  lines)  for  each  linear  polarization. 
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Dover, Gibraltar:  VV-Pol,n=0.026415,  c=0.36628,  n=459,  %Under30=14.8148,  %Over30=1 3.7255 


Dover, Gibraltar:  HH-Pol,|i=0.093396,  o=0.49449,  n=483,  %Under30=13.2505,  %Over30=1 6.9772 


AL/L  AL/L 


Figure  35  Histograms  of  AL/L  for  ships  at  incidence  angles  >  25°  showing  the  30%  error 
positions  (red  lines)  for  each  linear  polarization. 
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VUSAR  ship  signature  length  [m]  VUSAR  ship  signature  length  [m] 


Dover, Gibraltar:  HH-Pol,  n=654 


Dover, Gibraltar:  VV-Pol,  n=633 


Dover, Gibraltar:  HV-Pol,  n=373  Dover, Gibraltar:  VH-Pol,  n=392 


Figure  36:  Scatterplots  of  VUSAR  ship  length  and  ISR  Database  validated  ship  length  for  each 
linear  polarization.  The  blue  points  indicate  ships  with  TCR  >  10  dB  while  the  red  points 

indicate  ships  with  TCR  <10  dB. 
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Dover, Gibraltar:  HH-Pol,(i=-0.3508,  o=0.45585,  n=175,  %Under30=63.4286,  %Over30=9.1429 


AL/L 


AL/L 


Dover, Gibraltar:  HV-Pol,  (1=0.47171, 0=0.94195,  n=29,  %Under30=17.2414,  %Over30=41 .3793 


AL/L 


Dover, Gibraltar:  VH-Pol, (1=0.048309,  o=0.50087,  n=21,  %Under30=28.5714,  %Over30=33.3333 


AL/L 


Figure  37:  Histograms  of  AL/L  for  ships  with  TCR  <  10  dB  showing  the  30%  error  positions  (red 

lines)  for  each  linear  polarization. 
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Dover, Gibraltar:  HH-Pol, (1=0.085731,  0=0.49266,  n=479,  %Under30=1 3.7787,  %Over30=1 6.2839 


Dover, Gibraltar:  VV-Pol, (1=0.00026094,  0=0.35516,  n=476,  %Under30=16.1765,  %Over30=1 0.5042 


Figure  38:  Histograms  of  AL/L  for  ships  with  TCR  >  10  dB  showing  the  30%  error  positions  (red 

lines)  for  each  linear  polarization. 
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Dover, Gibraltar:  HH-Pol,  Error  Rate  =  40.2141%,  n=654 


Dover, Gibraltar:  HV-Pol,  Error  Rate  =  30.2949%,  n=373 


50  100  150  200  250  300  350  400  450 

Validated  ship  length  [m] 


Figure  39:  Scatterplots  ofVUSAR  ship  length  and  ISR  Database  validated  ship  length  for  each 
linear  polarization.  The  blue  points  indicate  range  travelling  ships  while  the  red  points  indicate 

azimuth  travelling  ships. 
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AL/L 


Dover, Gibraltar:  HV-Pol,|i=0.040628,  0=0.33638,  n=234,  %Under30=7.6923,  %Over30= 10.2564 


-0.5  0  0.5 


.5  2 


Dover, Gibraltar:  VH-Pol,n=0.087599,  o=0.37036,  n=318,  %Under30=6.6038,  %Over30= 14.7799 


Figure  40:  Histograms  of  AL/L  for  range  travelling  ships  showing  the  30%  error  positions  (red 

lines)  for  each  linear  polarization. 
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Dover, Gibraltar:  HV-Pol,  (1=0.59425,  a=0.94558,  n=139,  %Under30=4.31 65,  %Over30=48.9209 
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AL/L 


Dover, Gibraltar:  VH-Pol,(i=0.6874,  o=0.971 ,  n=74,  %Under30=4.0541 ,  %Over30=51 .3514 


Figure  41:  Histograms  of  AL/L  for  azimuth  travelling  ships  showing  the  30%  error  positions  (red 

lines)  for  each  linear  polarization. 
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7  Ship  wakes 


It  is  well  known  that  ship  wakes  appear  in  SAR  ocean  images,  and  that  ship  wakes  are  most 
visible  under  moderate  wind  conditions  when  the  clutter-to-noise  ratio  is  relatively  high.  A  new 
observation  from  this  data  set  is  that  there  is  a  signature  related  to  ship  wakes  in  the  cross¬ 
polarization  data. 

Figure  42  and  Figure  43  show  the  nature  of  the  observed  ship  wakes  in  both  co-polarization  and 
cross-polarization  data.  This  IS3  image  was  acquired  on  27  November  2006  with  HH  and  HV 
channels.  In  the  figures,  the  blue  line  is  the  projected  ship  track  as  calculated  from  the  A1S  data; 
the  intersection  of  the  blue  and  yellow  lines  represents  the  AlS-projected  ship  position;  the 
opposite  end  of  the  yellow  line  at  or  near  the  ship  target  signature  position  represents  the  AIS- 
predicted  ship  position,  which  accounts  for  azimuth  shifting  due  to  the  radial  component  of  ship 
velocity.  The  blue  lines  lie  along  several  wakes  that  appear  as  linear  features  in  the  HH  data.  The 
‘zooms’  illustrate  that  the  ship  signature  is  deterministically  shifted  in  azimuth  with  respect  to  the 
ship  wake  by  an  amount  that  is  predicted  by  the  length  of  the  yellow  line.  Of  particular  interest  is 
the  ghost  ship  signature  that  is  present  for  some  ships  in  the  HV  data.  The  ghost  signature  is 
thought  to  be  related  to  the  ship  wake  since  the  signature  is  not  subject  to  the  azimuth  shift 
associated  with  the  radial  component  of  the  ship  velocity.  Due  to  low  clutter  levels,  wakes  are 
not  expected  to  be  visible  in  cross-polarized  data.  However,  it  is  possible  that  the  ghost  signature 
could  be  related  to  scattering  off  the  bow  wave  of  the  ship. 

All  of  the  Envisat  ASAR  images  acquired  were  scanned  for  the  occurrence  of  ship  wake 
signatures  in  each  available  polarization  channel.  It  was  found  that  linear  ship  wake  signatures 
were  not  visible  in  approximately  50%  of  the  co-polarization  images  that  were  acquired.  Of 
course,  the  co-polarization  ship  wakes  were  not  necessarily  visible  for  every  ship  in  the  image 
that  included  ship  wakes.  The  co-polarization  ship  wakes  were  more  visible  for  smaller  incidence 
angles  when  the  clutter  (i.e.,  wind  speed)  wasn’t  too  high.  On  the  other  hand,  some  cross¬ 
polarization  wakes  were  visible  in  all  of  the  images  that  were  acquired.  Again,  the  cross¬ 
polarization  ship  wakes  were  not  necessarily  visible  for  every  ship  in  a  particular  image. 


46 


DRDC  Ottawa  TM  2008-005 


Figure  42:  Ship  wakes  observed  in  the  IS3  Dover  image  acquired  on  27  November  2006. 


Figure  43:  Zooms  of  ships  A,  B,  and  C  of  Figure  42  in  both  HH  (left)  and  HV  (right) 

polarizations. 
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8  Conclusions 


We  have  successfully  used  AISLive  data  to  validate  ship  signatures  in  Envisat  ASAR  AP  mode 
imagery.  The  data  sets  were  acquired  over  high  density  shipping  in  the  Strait  of  Dover  and  the 
Strait  of  Gibraltar.  This  work  extends  previous  studies  that  focussed  on  R-l  Fine  and  SCNB 
modes. 

In  total,  35  dual  polarization  images  were  acquired.  These  images  included  672  validated  ships 
for  HH  polarization,  658  for  VV,  387  for  HV,  and  396  for  VH.  The  AISLive  validation  data  were 
verified  through  the  1SR  database,  using  the  1MO  number  as  a  cross  reference.  The  mean 
distance  between  the  AlS-predicted  ship  position  and  the  actual  ship  signature  position  was 
around  150  m,  which  is  roughly  the  size  of  the  average  ship  in  the  validation  database. 

This  data  set  has  provided  new  insight  to  ship  detectability  in  SAR  ocean  imagery.  For  example, 
at  co-polarization  and  for  smaller  incidence  angles,  ships  may  not  be  visible  at  all.  On  the  other 
hand,  the  ship  is  usually  clearly  visible  in  the  corresponding  cross-polarization  channel. 

For  each  validated  ship,  several  metrics  were  computed  using  the  VUSAR  tool.  These  included 
the  clutter  statistics,  the  total  RCS  of  each  ship,  and  the  ship  signature  length.  These  metrics  were 
evaluated  in  terms  of  the  polarization,  the  local  incidence  angle,  the  validated  ship  length,  the 
target-to-clutter  ratio,  and  the  ship  aspect  angle.  The  following  is  apparent: 

•  For  co-polarization,  the  ocean  clutter  decreases  with  increasing  incidence  angle.  For  cross¬ 
polarization,  the  ocean  clutter  rarely  exceeds  the  instrument  noise  floor.  Therefore,  the 
cross-polarization  clutter  is  at  the  Envisat  ASAR  noise  floor,  so  any  contrast  metrics  would 
be  relative  to  the  instrument  noise  floor  rather  than  the  ocean  clutter. 

•  For  all  linear  polarizations  considered,  the  clutter  or  noise  may  be  characterized  as  K- 
distributed.  The  cross-polarization  clutter  tended  to  be  more  Gaussian  than  the  co¬ 
polarization  clutter. 

•  Several  ship  target  RCS  metrics  were  considered.  Straightforward  model  fits  provided  a 
relationship  between  the  total  RCS  and  the  ship  length.  The  HH  and  VV  total  RCS  are 
comparable  and  similar  to  previous  measurements  made  with  R-l  HH  data  and  with  EC  CV- 
580  C-band  polarimetric  SAR  data.  The  cross-polarization  total  RCS  is  about  10  dB  smaller 
than  the  co-polarization  total  RCS. 

•  There  is  significant  variability  in  the  RCS  and  contrast  metrics  considered.  Focussing  on  the 
segmented  region  RCS-to-clutter  ratio,  it  is  apparent  that,  for  co-polarization  the  contrast 
increases  with  increasing  incidence  angle.  This  is  because  the  ocean  clutter  decreases  with 
increasing  incidence  angle  for  co-polarization  data.  On  the  other  hand,  for  cross¬ 
polarization  the  contrast  is  more-or-less  independent  of  incidence  angle.  This  is  because  the 
noise  floor  against  which  the  contrast  measure  is  made  is  essentially  independent  of 
incidence  angle. 

•  The  ship  target  contrast  is  generally  larger  for  larger  ships. 
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•  The  ship  target  contrast  is  best  for  cross-polarization  at  smaller  incidence  angles,  and  is  best 
for  co-polarization  at  larger  incidence  angles.  The  break  point  between  the  two  cases  does 
depend  upon  the  ship  length,  but  occurs  at  about  33°  incidence  angle. 

•  The  total  RCS  variability  was  successfully  modelled  by  fitting  the  normalized  total  RCS  to 
several  well  known  PDFs.  The  RCS  variability  is  smaller  for  cross-polarization  ship 
signatures  than  for  co-polarization  ship  signatures. 

•  The  estimated  ship  signature  length  agreed  with  the  validated  ship  length  for  60. 1  %  of  the 
co-polarization  ships  and  71.5%  of  the  cross-polarization  ships.  For  co-polarization,  the 
error  rate  is  higher  for  smaller  incidence  angles  (64.7%)  and  errors  tend  to  arise  from 
underestimates  of  the  ship  length.  Similarly,  the  error  rate  is  higher  for  smaller  TCRs 
(65.2%).  These  trends  are  attributed  to  higher  clutter  and  lower  target  contrasts  for  smaller 
incidence  angles  for  co-polarization.  For  cross-polarization,  the  error  rate  is  also  higher  for 
small  incidence  angles  (35.4%),  but  errors  tend  to  arise  from  overestimates  of  ship  length. 
Similarly,  the  error  rate  is  higher  for  small  TCRs  (60.3%),  however,  there  were  very  few 
ships  that  had  a  contrast  smaller  than  1 0  dB  in  the  cross-polarization  data.  These  trends  are 
attributed  to  the  generally  higher  target  contrast  providing  more  visibility  for  sidelobes  or 
other  imaging  artefacts  at  cross-polarization.  For  all  polarizations,  there  tends  to  be  more 
overestimates  in  ship  length  for  azimuth  travelling  ships. 

•  Ship  wake  signatures  appear  in  the  cross-polarization  channel  as  a  ship-like  signature  that 
has  not  been  subject  to  azimuth  shifting,  as  occurs  for  the  ship  signature  itself.  These 
signatures  appear  relatively  frequently  and  may  be  due  to  scattering  from  waves  generated 
off  the  bow  of  the  ship. 

This  is  a  very  rich  data  set  that  has  yielded  new  observations  and  insights  to  ship  detectability  by 
SAR.  These  results  could  have  bearing  on  the  design  of  future  SAR  modes  or  SAR  missions,  and 
could  be  used  to  improve  the  performance  of  ship  detection  software. 

The  use  of  AIS  data  for  ship  signature  validation  is  now  a  well-proven  methodology  that  should 
be  used  to  develop  ship  detection  models  for  other  missions  such  as  TerraS AR-X  and  Cosmos- 
Skymed.  Focussing  on  high  density  shipping  regions  would  provide  a  large  set  of  validated  ship 
targets  for  each  image  acquired. 
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List  of  symbols/abbreviations/acronyms/initialisms 


A 

Ascending 

AIS 

Automatic  Identification  System 

AO 

Announcement  of  Opportunity 

AP 

Alternating  Polarization 

ASAR 

Advanced  SAR  (on  Envisat) 

CDF 

Cumulative  Distribution  Function 

D 

Descending 

dB 

Decibel 

DRDC 

Defence  R&D  Canada 

EC 

Environment  Canada 

ECDF 

Empirical  Cumulative  Density  Function 

EOLI 

Earthnet  OnLine  Interactive 

ESA 

European  Space  Agency 

GUI 

Graphical  User  Interface 

HH 

Horizontal  transmit,  Horizontal  receive  polarization 

HV 

Horizontal  transmit,  Vertical  receive  polarization 

EO 

Earth  Observation 

IMO 

International  Maritime  Organization 

IS# 

Image  mode  #  (ASAR) 

ISR 

Internet  Ships  Registry 

K-S 

Kolmogorov-Smimov  (goodness-of-fit  test) 

MSVR 

Mean- Squared  to  Variance  Ratio 

PDF 

Probability  Density  Function 

R-l 

RADARSAT-1 

R&D 

Research  &  Development 

RCS 

Radar  Cross  Section 

SAR 

Synthetic  Aperture  Radar 

SCNB 

ScanS AR  Narrow  B 

TCR 

Target  to  Clutter  Ratio 

VH 

Vertical  transmit,  Horizontal  receive  polarization 

VV 

Vertical  transmit,  Vertical  receive  polarization 
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Ships  appear  in  synthetic  aperture  radar  (SAR)  ocean  imagery  as  bright  targets  against  the  ocean 
clutter  background.  Envisat  Advanced  SAR  (ASAR)  Alternating  Polarization  (AP)  mode  data 
acquisitions  over  the  Strait  of  Dover  and  the  Strait  of  Gibraltar,  two  high-density  shipping 
regions  with  Automatic  Identification  System  (A1S)  coverage  via  AISLive,  have  provided  a 
large  database  of  validated  ship  signatures.  For  each  validated  ship  signature,  several  metrics 
were  computed  including  the  clutter  statistics,  the  total  radar  cross  section  of  each  ship,  and  the 
ship  signature  length.  These  metrics  were  evaluated  in  terms  of  the  polarization,  the  local 
incidence  angle,  the  validated  ship  length,  the  target  to  clutter  ratio,  and  the  ship  aspect  angle. 
This  is  a  very  rich  data  set  that  has  yielded  new  observations  and  insights  to  ship  detectability 
by  SAR.  These  results  could  have  bearing  on  the  design  of  future  SAR  modes  or  SAR 
missions,  and  could  be  used  to  improve  the  performance  of  ship  detection  software. 
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